Interface mapping in two-dimensional random lattice models 
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We consider two disordered lattice models on the square lattice: on the medial lattice the random 
field Ising model at T = and on the direct lattice the random bond Potts model in the large-g 
limit at its transition point. The interface properties of the two models are known to be related by 
a mapping which is valid in the continuum approximation. Here we consider finite random samples 
with the same form of disorder for both models and calculate the respective equilibrium states 
exactly by combinatorial optimization algorithms. We study the evolution of the interfaces with the 
strength of disorder and analyse and compare the interfaces of the two models in finite lattices. 
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I. INTRODUCTION 

The presence of different type of quenched disorder 
could influence the cooperative behaviour of many-body 
systems in different measure. In this respect dilution 
or random ferromagnetic bonds, which are often called 
as Tc-disorder, have a relatively weak effect. The fer- 
romagnetic order in the pure system, which is present 
for T < Tc, is not destroyed by a small amount of Tc- 
disorder. Dilution has a dramatic effect only in the vicin- 
ity of a phase transition point, \T — Tc\/Tc ^ 1. If this 
transition is of second order even its universality class 
can be changed for any small amount of disorder—. 

A disorder perturbation with a stronger effect is rep- 
resented by random fields^. As an example we consider 
here the random field Ising model (RFIM) in which the 
fields are taken from a distribution with zero mean and 
variance /ig. In dimensions, d > 2, the phase transition is 
governed by a zero-temperature fixed point, the ground 
state of the system being disordered for ho > ho^dd) > 0, 
and there is magnetic long-range-order for ho < ho^dd). 
In two-dimensions we have ft-o,c(2) = 0, thus the mag- 
netic order in the ground state is destroyed by any finite 
random field^i^. The ground state of the RFIM is con- 
veniently characterized by geometric (or Ising) clusters^, 
which are formed by adjacent spins being in the same 
orientation. In 2d these clusters are homogeneous only 
up to a finite characteristic length, £g, which is called 
the breaking-up length of geometrical clusters^"—. In a 
coarse-grained picture the ground state consists of ran- 
domly oriented domains of typical size, ig, and as a con- 
sequence there is no long-range-order in the system. For 
weak disorder, h^ <^ J , J ~ I being the nearest neigh- 
bour coupling constant, the breaking-up length is large, 
and for h^ ^ \t diverges as^: 



of bond disorder is also strong and the physical phenom- 
ena is analogous to that of the RFIM discussed before^. 
As a prototypical model having a first-order phase tran- 
sition we consider the g-state Potts modeliSi in which the 
number of states is sufficiently large, q > qdd). In the 
pure system having homogeneous couplings {J ~ Jo) at 
the first-order transition point, T = T^ the ordered and 
disordered phases coexist, the latent heat is finite and 
there is a finite jump of the order-parameter. Introducing 
random ferromagnetic couplings, such that the fluctuat- 

./ - Jo 
mg part, k = 



Jo 



igiho) r^ cM^ho^) , 



(1) 



with a disorder dependent constant, A = 0{1). 

As far as first-order transitions arc concerned the effect 



has zero mean and variance Kq the 

latent heat and the jump of the magnetization is gener- 
ally reduced with increasing kq. In dimensions d > 2 for 
weak disorder, Ko^c{d) > kq > 0, the transition stays first 
order, whereas for strong disorder, kq > Ko.c{d), the la- 
tent heat is vanishing and the transition softens to second 
orderi^iii. The critical exponents of this transition are in- 
dependent of the strength of disorder, but they generally 
depend on the value of qii^r— li^iii. In the borderline case, 
'*o = i'^o,cid), there is a tricritical point with a new class 
of tricritical exponents, which are expected to be inde- 
pendent of the value of q and the actual fixed point is at 
q -^ ooiS.. In 2d the limiting value is ko,c{2) = and the 
phase transition is softened to second order by any weak 
bond disorderi^iii. 



To characterize the state of the system we consider the 
high-temperature series expansioni^, which in the large-q 
limit is dominated by a single diagram. In the pure case 
this diagram is either the fully connected graph, which 
happens in the ordered phase, T < T^ or the fully discon- 
nected graph, which is the case in the disordered phase, 
T > Tc. At the phase transition point these two basic 
graphs coexist. In the disordered model the dominant 
diagram at T = Tc contains both connected and discon- 
nected parts. In 2d the optimal diagram is homogeneous 
only up to a characteristic finite length, £d, which is the 
breaking-up length of the optimal diagram. For small 
disorder the breaking-up length is large and in the limit 



Ko 



it is divcrgenti^: 



^d{K.n) --exp(BKo^) 



(2) 



which is in the same form, as for the RFIM in Eq.(IT|). 
In the thermodynamic hmit the optimal diagram is a 
random composition of the two basic graphs of typical 
size, £d, therefore there is no phase coexistence in 2d and 
the transition is of second order. 

As described above the RFIM and the random bond 
Potts model (RBPM) have analogous physical properties. 
This analogy was first noticed by Cardy and Jacobser^i^ 
(CJ) and they have obtained an exact mapping, which 
is valid in the limit of d — > 2 and q — > oo, when the in- 
terface Hamiltonian of the two problems have the same 
type of solid-on-solid (SOS) description. From this map- 
ping follows that the energy exponents of the tricritical 
RBPM are equivalent to the magnetization exponents 
of the critical RFIM. Furthermore, analyzing the rcnor- 
malization group (RG) flow it was conjectured that the 
above mapping stays valid for d > 2, in particular at 
d = 3. This conjecture has been numerically testedi^ for 
the 3d RBPM in the large-g limit. 

The C J- mapping has consequences at d = 2, as far 
as interface properties of the two models arc concerned 
at large scales. This has already been mentioned with 
respect to the same form of the breaking-up lengths in 
Eqs.([T]) and ^. In this context one may ask the question 
how this mapping works in finite scales? In this paper 
we are going to perform a systematic study and consider 
both models with exactly the same form of disorder. For 
this we put the Potts model on the square lattice, whereas 
the Ising model on the medial lattice. In this way cou- 
plings of the Potts model and fields of the Ising model 
are on the same position. We study random samples 
with fixed boundary conditions, which prefer the oppo- 
site orientation for the RFIM and the different phases 
(basic graphs) for the RBPM at the two halves of the 
lattice. Using combinatorial optimization algorithm^Sii^ 
we calculate exactly the ground state of the RFIM and 
the optimal diagram of the RBPM and study their evo- 
lution with increasing the strength of disorder. We also 
compare the diagrams obtained in the two problems and 
study their fitting. 

The structure of the rest of the paper is the follow- 
ing. The RFIM and the RBPM are defined in SecHIl 
together with the asymptotic mapping. Evolution of the 
diagrams with the strength of disorder are presented in 
Sec lIIII whereas a comparison of the diagrams of the two 
models are shown in Sec lIVI Our paper is closed by a 
discussion in SeclVl 



II. MODELS AND ASYMPTOTIC MAPPING 

We consider a square lattice the sites of which are at 
(x,2/)withl <x< Land-(L-l)/2< y< (L-l)/2and 
both X and y being integer. In the corresponding medial 
lattice the sites are placed to the middle of the links and 



have coordinates (u,w) as {x,y + 1/2) and (x + 1/2, y), 
see Fig HI 




i,t/- 1 



FIG. 1: The square lattice (black) with the Potts spins and 
its medial lattice (red) with the Ising spins. 



A. Random field Ising model 

Spins of the random field Ising model, au.v = ±1, sue 
put on the vertices of the medial lattice. The Hamilto- 
nian of the model is given by: 



'Hi = — 2_^ Ou.vOu 



+ l/2,-u±l/2 



^ /i(u, w)cr„^t, . (3) 



The external field, h(u, v), is a random number which is 
parametrized as: 



h{u,v) = hoe{u,v) 



(4) 



where the distribution of e(u, v) has zero mean and vari- 
ance unity. We apply fixed spin boundary conditions: 
the Ising spins at the upper part of the boundary (v > 0) 
are fixed to -f 1, and at the lower part of the boundary 
{v < 0) they are fixed to —1. This model is studied at 
T ~ 0. For a given realization of the disorder we calculate 
the ground state exactly by a combinatorial optimization 
algorithm. 



1. SOS approximation 

For weak disorder, such that ig{ho) > L, the ground 
state is well approximated by two clusters with -1-1 and 
— 1 spins, respectively, and having an interface in be- 
tween. Furthermore we use the solid-on solid (SOS) ap- 
proximation, when the height of the interface at position 
u is given by a unique function, C('")- In this case the 
position of the interface is obtained by minimizing the 
SOS energy functionali^: 



EiO 



-'^oE E - E h("'-) 



(5) 



iij<C(«) ">?("); 



with the constrain that the length of the interface, mea- 
sured in the medial lattice is constant. 



B. Random bond Potts model 

The random bond Potts model is put on the vertices 
of the original square lattice and defined by the Hamil- 
tonian: 



Hv 



E 



J{x,y + l/2)S{sx,y,Sx,y+i) (6) 



Here s^ 



1,2, 



J{X + 1/2, y)S{Sx,y, Sx+l,y)] ■ 

, g is a Potts spin variable at site 
{x,y) and the couplings, J{x,y + 1/2) = J{u,v) > 
and J{x + 1/2, y) = J{u,v) > 0, are independent and 
identically distributed random numbers. Introducing 
the reduced temperature, T' = Tlnq, and similarly, 
j3' = 1/T' = /^/Ing, the partition function in the ran- 
dom cluster representation is given hy^: 

^ = E 9'^"" n ['^'"^"'"^ - 1] (^) 

GCB (u,u)£G 

with K{u,v) = j3'J{u,v). Here the sum runs over all 
subset of bonds, G Q E and c{G) stands for the num- 
ber of connected components of G. In the following we 
consider the large-g limit, where q^^'^''"^ ^ 1, and the 
partition function can be written as 



GCE 



0(G) = c{G) 



E ^("'^) (8) 
{u,v)eG 

which is dominated by the largest term, (f>* — maxc 4>{G), 
so thati^i2^ 



r 



subleading terms. 



(9) 



For random couplings with < K{u,v) < 1 wc use 
the form, K{u,v) ~ K{\ + k.{u,v)) and the k{u,v) 
random numbers are taken from a symmetric distribu- 
tion: P(k) = P(— k), with a variance Kg, thus k{u,v) = 
Koe{u,v). For this type of distribution the phase tran- 
sition point of the system follows from self-duality and 
given by^i; 

Kc^l- (10) 

Thus at the critical point the random couplings are 
parametrized as: 

P'j{u, v) = K{u, v) = -[l + Kae{u, v)] (11) 

and the strength of disorder is measured by kq. In our 
problem we use such type of boundary conditions, that 
the boundary couplings at the upper part of the lattice 
{v > 0) are strong k+ = 1, thus promoting an ordered 
phase, whereas at the lower part of the lattice {v < 0) 
these are weak, k~ = 0, thus favouring the disordered 
phase. For a given realization of the disorder we calculate 
the ground state exactly by a combinatorial optimization 
algorithm^. 



1. SOS approximation 

In the following we consider the model at the critical 
point in the weak disorder limit, when idino) > L. In this 
case we use the approximation that the optimal diagram 
consists of a fully connected part, denoted by Gc-, and a 
fully disconnected part, denoted by Gc, having an inter- 
face in between. The interface is approximated with the 
solid-on-solid (SOS) model: having a unique height C(-u) 
at position v^^^^. The lattice contains N = L"^ points, 
from which n{Gc) and n{Gc) are in Gc and in Gc, re- 
spectively. Similarly, out of the E = 2L{L — 1) edges, 
there arc e{Gc) in the subgraph, Gc- The cost function 
in Eq.Q is given by: 

= N+l-n(Gc)+ Yl {'^ + i^Q'^iu,v))/2 (12) 
{u,v)eG 

= N+l-n{Gc)-eiGc)/2+ J^ «oe(w,w)/2. 

{u,v)eG 

Now we use the microcanonical condition: 
Ea/; bonds e("'^') = and arrive to the cost func- 
tion: 



40(c) = Ko^ E - E b("'^) 



(13) 



i)J<C(M) t-'>C(u}l 



for a fixed value of n{Gc) — e{Gc)/2. This latter quantity 
is uniquely determined by the length of the interface, 
which is measured in the square lattice. 

Comparing the two expressions in Eqs.(l5]) and (|13[) we 
see that in the SOS approximation the position of the 
interface is obtained by finding the extremal value of the 
same expression, however with somewhat different con- 
strains. In both cases the length of the interface is fixed, 
however in the RFIM this length is measured in the me- 
dial lattice, whereas for the RBPM in the original lattice. 
In the following we study both models numerically and 
see how well the SOS approximation and thus the above 
relation is valid in finite systems. 



III. EVOLUTION OF INTERFACES WITH THE 
STRENGTH OF DISORDER 

In the numerical study we use bimodal disorder, such 
that eiu, v) ~ ±1 with the same probability. 



A. Random field Ising model 

For RFIM the evolution of the ground state with the 
strength of disorder, h^, for a given realization of the 
random numbers, e{u,v), is illustrated in Figl2] In the 
limit ft-o ^^ 0"*" there is a unique SOS interface between 
two oppositely magnetized clusters. With increasing Hq 
at /iQ.i the ground state will change. It still consists of 



two large clusters, however the interface have some over- 
hangs, thus it is no-longer of SOS type. With further 
increased disorder, for ho > /10.2, no longer two clusters 
exist, but some islands are formed in the big clusters. 
Finally, for strong disorder, h > /io,3j the ground state is 
represented by several disjoint clusters and one can not 
identify an interface in the systerr^. 





b) 











d) 








FIG. 2: Evolution of the ground state of the RFIM with in- 
creasing strength of disorder in a 32 x 32 lattice. Up and 
down spins are represented by dark (red) and light (green) 
regions, a) ho = 0.1275: Two clusters with an SOS interface, 
b) ho = 0.2525: Two clusters with an interface with over- 
hangs, c) ho = 0.3275: Internal islands are formed within the 
big clusters, d) ho = 4.0: Several disjoint clusters are formed. 



We have studied the behavior of the first relevant dis- 
order scale, /iQ.i, at which the SOS picture is broken 
down. Having several samples we have calculated its 
mean value, which is then repeated for different finite 
sizes, L. ho^i{L) is found to go to zero with L as a 
power-low, ho^i{L) ^^ L~'^, which is illustrated in FiglS] 
The exponent is found to be: uj ~ 0.36(5). Similarly, the 
second disorder scale, /io,2, is found to follow the same 
type of decay with L, with the same exponent, w, but 
with a different prcfactor. 

Our numerical results arc in agreement with the exact 
result, that for any finite /iq > in the thermodynamic 
limit there is no long range order in the 2d RFIM. 



B. Random bond Potts model 

We have repeated the previous analysis for the RBPM. 
For a given random sample we have varied the strength 
of disorder, kq, and studied the form of the optimal di- 
agram. Having the same set of e(u, v) as for the RFIM 
in Figl5] the results are collected in FigH) In the limit 
kq ^ the interface is plane, with a minimal length, 
Ljnt = L (not shown in Fig|3|). At a limiting value, kq.Oj 
the interface becomes an SOS interface with Lint > L. 
Further increasing kq the length of the interface is in- 




FIG. 3: Average value of the critical disorder strength for the 
RFIM, /io,i, and that for the RBPM, ko.i, as a function of L 
in log-log scale. The straight line indicating the theoretical 
prediction has a slope — cj = —1/3. 



creasing, but it stays SOS type. Afterwards, at kq.i the 
interface will loose its SOS character: overhangs and bub- 
bles appear. Finally, for strong disorder the optimal di- 
agram is composed of disjoint subgraphs. We can thus 
conclude that the evolution of the interfaces in the RBPM 
and in the RFIM with the strength of disorder are quali- 
tatively very similar. There are, however, differences, for 
1^0 '^ 1^0, 0, which are due to the different orientations of 
the interface lines. 

We have also studied the size dependence of the aver- 
age value of Ko,i, which is shown in FiglJ] Here also a 
power-low decay is found with the same exponent, w, as 
for the RFIM. 

We can explain the value of the lo exponent, which is 
measured both for the RFIM and for the RBPM, in the 
following way. Let us concentrate now on the RFIM, 
for the RBPM similar reasoning works. At ft,o,i the 
first overhang in the ground state of the RFIM appears. 
The difference between the ground states for h < /io,i 
and h > ho^i has a narrow shape. Its typical width is 
w ~ 0(1), since the total length of the interface is in- 
creased by ALint = 0(1), and its height is given by L±, 
which is the size of the transverse fluctuations of the SOS 
interface and given by2^: L± ^ L^^^. Now we use the 
Imry-Ma argument^ and compare the loss of energy due 
to the increase of the interface: Eint ~ w, with that of 
the gain due to disorder fluctuations: Efi ^ ho^iy/'wL±. 
From the condition: Eint = Efi we obtain the estimate 
/io,i ^ L^^'^, thus the exponent is given by w = 1/3 in 
agreement with the numerical results. 



IV. COMPARISON OF THE STRUCTURE OF 
THE GRAPHS IN THE TWO MODELS 

The SOS mappings presented in Sec |lT] assure that for 
weak disorder the scaling behavior of the interfaces in 



d) 



FIG. 4: Evolution of the optimal set of the RBPM with in- 
creasing strength of disorder in a 32 x 32 lattice having the 
same set of e(it, i)) as for the RFIM in Fig[2] Connected 
graphs and isolated sites are represented by dark (red) and 
light (green) regions, a) kq — 0.25: Two clusters with an SOS 
interface, b) kq = 0.34: Two clusters with an interface with 
overhangs, c) kq = 0.5: Internal islands are formed within 
the big clusters, d) kq — 0.76: Several disjoint clusters are 
formed. 



the two models are asymptotically equivalent. Here we 
study the question what happens for finite systems and 
for not too weak disorder, which type of similarity oc- 
curs between the ground states of the two problems. For 
this study wc consider exactly the same random samples 
for the two models, which means that the set of random 
variables, e{u,v), are the same for each position, {u,v), 
just the strength of disorder, ho and kq, respectively, can 
be different. For a given sample we calculate the ground 
state of the RFIM, as well as the optimal set of the RBPM 
and compare them. In order to define a quantitative mea- 
sure of the difference between the two graphs we consider 
the medial lattice and assign to each site, («,«), a vari- 
able denoted by D{u,v). If in the ground state of the 
RFIM a{u, v) ~ +1 {(j{u, v) = —1) and at the same time 
in the optimal set of the RBPM the edge (u, v) is occupied 
(non-occupied), then D{u,v) = 0, otherwise D(u,v) = 1. 
In FiglS] we compare the ground state configurations in 
Fig 12] with the optimal sets in FigU) One can see in this 
figure that close to the boundaries D{u,v) is typically 
zero (Greyscale (gsc) 3 (red) and gsc 2 (green) points) 
and in the interface region we have sites with I?(u, v) = 1 
(gsc 4 (blue) and gsc 1 (yellow) points). The difference 
between the two graphs is defined by the fraction of non- 
coherent sites: 



S{ho,Ko) 






2L{L-1) 



(14) 



what we call as discrepancy. Here 2L(L— 1) is the number 
of sites in the medial lattice with the given boundary 
condition. 

In the actual calculation we have fixed the value of 



the RBPM disorder parameter, kq, and calculated the 
discrepancy for different values of /ig. In this way we 
have measured the minimum value of the discrepancy, 
Smini^o), and the corresponding value of the RFIM 
parameter—, ho. The average value of the minimal dis- 
crepancy, [i5mm]av('^o) vcrsus Kq is plotted in Figini for 
different finite systems from L = 16 to L = 128. 



FIG. 5: Comparison of the ground state of the RFIM with 
the optimal set of the RBPM for the same pair of disorder 
parameters {ho and kq) as given in the panels of Figs[2]and 
m respectively. Greyscale (gsc) 3 (red) points represent: up 
spins in the RFIM and the connected graph in the RBPM; 
gsc 2 (green) points: down spins and isolated sites; gsc 4 
(blue) points: down spins and connected graph; gsc 1 (yellow) 
points: up spins and isolated sites. The gsc 4 (blue) and 
gsc 1 (yellow) points are non-coherent sites. The values of 
the discrepancies are S = 24/1984, 135/1984, 187/1984 and 
559/1984, for panels a), b), c) and d), respectively. 



For large kq the curves have a plateau, then with de- 
creasing Kq they start to decrease, pass over a minimum 
and afterwards increase for small kq. With increasing 
size the value at the minimum [<5mm]™"(i) is decreasing 
and in the largc-i limit the minimum is expected to be 
shifted at the origin. Indeed in the inset of FigjS]wc have 
plotted [(5mm]™"(i) as a function of L and in a log- log 
scale an asymptotically linear dependence is found, thus 
[<5™m]™"(i) ~ L-'^'. Here the exponent is w' = 0.39(7), 
which is compatible with w' = 1/3, the value of which fol- 
lows from the following simple argument. For weak disor- 
der the interface in the RBPM is approximately fiat (see 
the reasoning in Scc lIIip whereas in the RFIM it is rough 
having a transverse fluctuation of size L± ^ L"^/^. In the 
interface region there are L x L±^ sites, consequently the 
average discrepancy behaves as ^ L~^/^, in agreement 
with the numerical results. We can thus conclude that 
in the large-L limit the behavior of the discrepancy with 
the strength of disorder is the following. It starts from 
zero at Kq =0, has an approximately quadratic variation 
for small disorder and approaches a plateau for stronger 
disorder. This behavior is in complete agreement with 
the SOS mapping in SccHH. 



We have also checked the actual form of the graphs in 
the two problems for strong disorder. As illustrated in 
Fig l5] the two graphs generally have quite similar struc- 
ture, as far as the subgraphs, the topology of the islands, 
etc. are concerned. 




FIG. 6: Average value of the minimal discrepancy as a func- 
tion of the disorder parameter kq for different finite systems. 
In the inset the minimal value of [(5mm]av('«o) is plotted as 
a function of L in a log-log scale. Here the straight line has 
a slope of oj' = 0.39, to be compared with the theoretical 
prediction: f/3. 



V. DISCUSSION 

In this paper we have considered two problems of lat- 
tice statistics in which disorder has a strong and some- 
what analogous effect as far as the cooperative behavior 
of the models are considered. In both problems geomet- 
rical interpretation of the state of the system is used. 
In the first problem, which is the two-dimensional Ising 
model we use clusters of parallel spins (geometrical clus- 
ters) to characterize the ground state. Here random fields 
destroy the ferromagnetic order at T = and the geo- 
metrical clusters have homogeneous parts of finite ex- 
tent, Ig^ho). In the second problem, which is the two- 
dimensional Potts model with a large number of states 
the optimal diagram of the high-temperature series ex- 
pansion is used as a geometrical interpretation. Here 
we consider the first-order transition point of the pure 
system, which transition is softened to a continuous one 



due to bond disorder. As a consequence at the transi- 
tion point there is no phase-coexistence and the homoge- 
neous parts of the optimal diagram have a finite extent 
oiidiKo). 

In this paper we have considered finite samples of linear 
length, L, by varying the strength of disorder, Hq and 
kq in the two problems, respectively. Having fixed spin 
boundary conditions for £g{ho) > L {£d{i^o) > L) there 
are basically two phases (two elementary diagrams) in the 
ground state (optimal diagram) of the RFIM (RBPM), 
which are separated by an interface. For large scales and 
in the SOS approximation the interface Hamiltonians of 
the two problems have similar forms. 

Here we have studied the validity of the interface map- 
ping in finite systems, by putting the Potts model on the 
square lattice and the Ising model on the medial lattice. 
In this way the bonds of the Potts model and the fields of 
the Ising model have the same location. Using the same 
disordered samples we have studied the evolution of the 
interfaces in the two models as the strength of disorder is 
gradually increased. We have seen a similar trend of the 
evolution in the two models, however, in finite lattices 
we have observed also differences. These differences are 
basically due to the fact, that the interfaces have differ- 
ent orientations in the two models: in the RBPM it is in 
the (1,0) direction, whereas in the RFIM it has an (1,1) 
orientation. In the SOS approximation the interfaces are 
obtained by minimizing the same cost function, however 
with different constrains in the two problems. 

We have also compared the cluster structure of the 
two problems and have obtained the following conclusion. 
The relative difference between the two geometrical ob- 
jects as quantified by the discrepancy in Ea.(|14p tends to 
zero if i) the size of the system goes to infinity and ii) the 
strength of disorder goes to zero. The finite-size correc- 
tions are found to be in power law form, which vanish as 
L~" with uj = 1/3. Consequently relatively large finite 
samples are needed to see the asymptotic behavior. 
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